Behavior of tumors under nonstationary theraphy 
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We present a model for the interaction dynamics of lymphocytes-tumor cells population. This 
model reproduces all known states for the tumor. Futherly, we develop it taking into account 
periodical immunotheraphy treatment with cytokines alone. A detailed analysis for the evolution 
of tumor cells as a function of frecuency and theraphy burden applied for the periodical treatment 
is carried out. Certain threshold values for the frecuency and applied doses are derived from this 
analysis. So it seems possible to control and reduce the growth of the tumor. Also, constant values 
for cytokines doses seems to be a succesful treatment. 
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I. INTRODUCTION 

Cancer is one of the leading research areas, since this 
desease is a main cause of death. Surgery and chemother- 
apy are unsuccessful in many cases. Today the principal 
efforts are addressed to search new treatment strategies 
e.g. in immunotherapy (see Refs. Jl]]^] and references 
therein). In this case we refer to the use of cytokynes 
to stimulate the immune system. This is a protein hor- 
mone produced by activated lymphocytes which mediate 
both natural and specific immunity. The use of cytokines 
alone to boost the immune system represent one of the 
methods more commonly used in immunotherapy. The 
temporal evolution for this treatment comprises differ- 
ent steps: by supplying a starting dose of cytokines the 
rate of lymphocytes begins to increase due to the im- 
munological reaction Q reaching a maximum value; af- 
terwards the lymphocytes begin to decay because of the 
decrease of cytokines concentration inside the body until 
it reaches normal values. In the following we shall call 
this time the activation period of the immune system. 
This proccess repeats again between any two succesive 
dose suplies. We consider these injections separated in 
time by the dosage period. Otherwise, it would provoke 
an overdose or failure in the treatment for shorter and 
longer time, respectively. 

This work is devoted to understand the temporal tumor 
behavior when a periodical immunotherapy treatment is 
provided. Besides, we want to explore the set of parame- 
ter values to reproduce with our model the same features 
presented in | (like short tern oscillation in tumor size 
and long term tumor relapse). 

For this reason we reformulate the predator-prey model 
including new terms in the model which give account of 
tumor agressiveness, the diffusion of lymphocytes and the 
effect caused by cytokines on the tumor (Section 2). We 
analyse the model from the mathematical and biological 
points of view (Sections 3 and 4). We study in detail 



tumor evolution for different treatment regimes (Section 
5) and, in the last section, we discuss our results and 
suggest some options for improving the model. 



II. THE MODEL 

An extensive review of all models of tumor-immune 
system dynamics p|-|l2| as far as we know was done. We 
agree with the idea that such dynamics is determined 
by a competence of interacting species resembling the 
predator-prey model. The main difference between our 
model and the clasical deterministic model (e.g. Bell 
1973) B is the inclusion of new terms taking into ac- 
count, (1) the death of lymphocytes due to the increase 
of malignant cells population, (2) the flux of lymphocytes 
towards the place of local interaction and (3) the effect 
produced by the application of cytokine doses. 
Let X and Y denote respectively the number of malig- 
nant and lymphocyte cells. The rate of malignant cells 

(if) is S iven b y : 



dX 
~~dt 



= aX - bXY 



(1) 



We assume a growth rate proportional to X and a de- 
crease rate proportional to the frecuency of interaction 
with lymphocytes. Coefficients are a and 6, respectively, 
where a is tissue dependent. 

On the other hand, the growth rate of lymphocytes 



is described by: 



kX + u 



(2) 



It is proportional to the interaction with malignant cells 
and also to the flux per unit time of lymphocytes to the 
place of interaction. These effects are represented by the 
first and fourth terms in the right-hand side of equation 
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g. This last term characterizes the difussion process of 
lymphocytes that takes place in the surroundings of the 
tumor assuming a constant lymphocytes flux . On the 
other hand, the decrease rate depends on two factors: 
natural death and growth of malignant cells related to 
the effective area of tumor interacting directly with the 
lymphocytes. These are given by the second and third 
terms of this same equation where / and k are their re- 
spective coefficients of proportionality. 

In order to introduce the effects produced by the treat- 
ment with cytokines in the proccess of activation of the 
immune system, we add a periodical function that mim- 
ics the periodical dosage. As a first approximation, we 
propose the function Fcos 2 cut, where oj is the frecuency 
of the periodical behavior for cytokines inside the body. 
The modified Eq ^| will be given by the expression 



dY 

— = dXY - fY- kX 
dt 



F cos tut 



(3) 



Taking into account Eq.|l], we get the following system 
of differential equations. 



dX 

~dt 

dY 

~~dt 



aX - bXY 
dXY - fY + Q 



Q = - KX + u + Fcos 2 (ojt) 



(4) 

(5) 
(6) 



where Q is a function that picking up all the news con- 
tributions respecting the standard predator-prey model, 
with X(0) = Xq and Y(0) = Yq as initial conditions. 

From the analysis for Q = the system of Eq^, 
Eq.[|and Eq.|| reduces to the predator-prey model, which 
has as equilibrium points a saddle point at the phase 
portrait origin and a center in the first quadrant of the 
phase diagram p|. This center reflects the oscillating 
behavior of the competence between both predator and 
prey, which is characterized by an oscillation frecuency 
(number of cycles per unit time around the center of the 
phase diagram). Taking the reverse of this frecuency as 

a characteristic time to = —j= and rescaling the equa- 

v a f 

tions of the system of Eq.[|, EqJ|and EqJ| by means of 
the following scaling parameters: 

t = t r 
X = X'x 
Y = Y'y 

the system given by the Eq^, Eq.||and EqJ| become: 



VI 

— x= axX' - bxyX'Y' 
to 

~to~ 



y = (dxX'- f)yY' + Q 
Q = Fcos 2 (u t T) + u-K X'x 



(7) 

(8) 
(9) 



sustituting the scaling parameters by the following 
values: 



to = 
X' = 
Y' = 



1 



^7 



d 

/57 



we get the following rescaled equations 
dx 

dy 



a x — xy 
1 



xy 



y- kx + a + V cos (/3 r) 



(10) 
(11) 



with x(0) — xo and y(0) — yo as initial conditions. 

where V=%, k = f d VaJ, 

a = ./%, a = 4 and [3 = -£= 

Analysing in detail the prior expressions we can interpret 
these parameters as follows: 

From V it is deduced that its value depends on the 
net value for doses (F) and on the action of the immune 
system on malignant cells (b). Hence it would represent 
an effective value of all doses employed in the activation 
of lymphocytes, namely, those necesary for activating the 
lymphocytes that take action directly against the tumor 
cells. 

On the one hand, k is directly proportional to K which 
accounts for the negative effects exerted on the popula- 
tion of lymphocytes due to the size of the tumor and it 
is inversely proportional to the recognition (d) and at- 
tack (b) frequencies of the immune systen to malignant 
cells. So, we infer that the inverse value of this parameter 
(1/fc) gives, "in some sense", the control exerted by the 
immune system over the aggressiveness of the tumor due 
to its size. 

From the similarity of the expressions for a and V, taking 
into account the remarks done for V we interpret that 
a which depends on the action of the immune system 
on malignant cells (b) represent of all lymphocytes flux 
u those effective value involved in the attack on tumor 
cells. Finally, a and (3 are directly related to the prolifer- 
ation of malignant cells and the frecuency of treatment, 
respectively. 



III. DYNAMICS AND STABILITY ANALYSIS 
WITHOUT TREATMENT 

The system of Eq.[l0| and Eq.[ll] corresponding to V = 
0( no treatment) is an autonomous system given by the 
following equations: 



xy 



dx 

d y 1 fr m 

— = xy y — kx + a 

dr a 



(12) 
(13) 
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with x(0) = xq and y(0 ) — j/o as initial conditions. 
Susbstituting the Eq.[12] in the Eq.[l3| we get the differen- 
tial equation 

d x , 1 1 dx , dx , , , o , , , 

1-9 +( x i-)~r = (k- a)x 2 +ax 14) 

dT z a x ar ar 

with x(0) — xq, x(0) — vo as initial conditions. 
This equation is similar to that describing the motion 
of a particle in a force field 113] , whose potential is: 



U(x) — —-(k — a)x 



1 



(15) 



This potential has two extremes given by 

and — - a ~ x 



x\ 



X-2 



These extreme points depend on a, k and a as: 
X\ = minimum 



(7 > 1 



and for 



a < 1 



- > 1, 



- < 1, 

>1, 



< 1, 



x 2 > 
x 2 < 



maximum 



maximum 



maximum 



x 2 > 



minimum 



x 2 < minimum 
We only consider motion for x > 0, the suitable potential 
fields describing the motion of the particle are depicted in 
Fig. la and Fig. lb. In Fig. la , the maximum represents 
an unstable point for the particle motion, contrary in Fig. 
lb the particle oscillates around the minimum. 
The analysis of fixed points in the phase space for Eq.|l2] 
and Eqjl3| shows two steady-states. A fixed point is Lq = 
(0, aa) with associated eigenvalues 



a 2 (l -a) -I 
2a 



± 



« 2 (1-ct) + 1 



2a 



(16) 



For a < 1, Lq is a saddle point while for a > 1 is a stable 
node. 

The other fixed point is L% = (jjjE% j a ) with associated 
eigenvalues given by: 



A± 



aa 



2a(a — k) 



± 



2a(a — k) 



-(1-a) (17) 



The real part of this eigenvalue is zero for — = a with 
u <l. 

When condition (a — k) 2 — — > — 1 is fulfilled, we get two 
values for a c solutions of 

a 2 a 2 - 2a[k - 2a(a - k) 2 ]a c + k 2 - Aa 2 (a - k) 2 = 
given by 



a c = - - 2a 2 ( 1 - - 
a \ a 



±2a|l--U/a 2 

ay V a 



k 

- + 1 

a 



(18) 



defining the region of complex eigenvalues and focus- like 
behavior. 

The analysis of these eigenvalues provides a rich dy- 
namics. For the case — < 1 we have different situations. 
The states and its stability for the second fixed point are 
depicted schematically in Fig. 2. 

If a < | (ReX± > 0), we have an unstable focus or node 
depending on the parameter value a c relative to those 
given by Eq.|l8|. 

On the contrary, stable behavior (focus or node) appears 
when | < a < 1 (ReX± < 0). 

Now, if 1 < a the fixed point corresponds to a negative 
population of malignant cells, with no physical meaning. 
For the case £ > 1 and a < 1 < £ (Re\± < 0) the 
critical point moves to the second quadrant of the phase 
diagram, being discarded as before. 

For values of a in the ranges 1 < a < ^ and 1 < ^ < a we 
get a saddle point (A+ > and A_ < 0), whose separatrix 
splits the phase portrait into stable and unstable zones 
as can be seen in Fig. 5. In all the cases the dynamics 
in the phase diagram is represented by a homeomorfism 
fL4f between two fixed points. 



IV. BIOLOGICAL SIGNIFICANCE 

So far, we have presented a detailed analitical study 
of the linear stability of our model when V is set equal to 
zero. The interpretation of this preliminary results will 
give us the esential features of the system. 
Let us start with the case — < 1. 

a 

For a < — the system evolves towards a state of uncon- 
trolable tumor growth(see Fig. 3a) This case can be inter- 
preted as a recurrence like behavior J^|,[h| very similar 
to q-switching oscillations observed in physical phenom- 
ena as, for example, in lasers. On the contrary, when 
— < a < 1, our system evolves towards a controlable 
mass of malignant cells in a damped oscillating way (Fig. 
3b). This state is considered by some authors as a dor- 
mant state 

However in both cases, there exist populations of malig- 
nant cells that grow towards a state in which immunolog- 
ical activity has been suppressed. In the first case this 
happens for any initial conditions, while in the second 
it only happens for an initially weak inmunological re- 
sponse (Fig. 4a and 4b). Let us now analyse the reverse 
condition — > 1. 

a 

In this case there are two possible ranges for a: 
1 < - < a and 1 < a < - 

a a 

In both situations we are in the presence of a saddle 
point which means that for populations of cancer cells be- 
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low the horizontal separatrix the dynamics is irreversible: 
this curve represents the critical amount of malignant 
cells for a fixed population of lymphocytes. 
This situation is similar to the case analised before for a 
weak immnune system as an initial condition: the popu- 
lation of malignant cells grows towards a value such that 
the immunological response is reduced to zero (Fig. 5). 
This would represent a state where illness is not the cause 
of death but leaves the body unprotected against other 
diseases. 

However for malignant cells above the horizontal separa- 
trix it is possible to observe regression of tumor as has 
been reported in clinical experiments (see Ref. jl5| and 
references therein) 



V. STABILITY ANALYSIS WITH TREATMENT 
AND BIOLOGICAL IMPLICATIONS 



The system represented by Eq.|l0| and Eq.ll with 
F / (cytokines doses amplitude) can be analised as 
an autonomous system fl7|| . The procedure consists in 
substituting the oscillating function cos (3t of the driven 
term F cos 2 j3t in the Eq.Qh by a new variable u, which 
is a solution of the second order differential equation of 
a linear oscillator 



^ +/3 M = 



(19) 



(where it(0) = 1, u(0) = 0) which can be written as two 
linear coupled differential equations 



dz 
du 



-(3 2 u 



(20) 
(21) 



(with it(0} = 1, z(0) = 0) 
Then Eq.[n^ and Eqjll] become: 



xy 



xy 



-p 2 u 



— ka 



(22) 
(23) 
(24) 
(25) 



with x(0) = Xq, y(0) = ijq, u(0) = 1, z(0) = as initial 
conditions. 

This system presents the critical points Lq — (0, aa, 0, 0) 

and L\ = f^E^i a, 0, 0^j whose projection in the y-x 

plane coincides with those critical points of the unper- 
turbed system (Eq. [12|a nd Eq. |13|) with the same eigenval- 
ues given by the Eq.|l6| and Eqjl7|plus the new conjugate 
pair A± = ± i/3. 



In this case we are in presence of a center manifold 
where solutions can be expanding or contracting, i.e., the 
asymptotic stability analysis carried out before loses its 
validity, needing more complex developments. 
In order to avoid such complexity and gain a better com- 
prehension, we may consider this system like a couple of 
one linear (Eq.[l9|) and one nonlinear (Eq.|l4|) oscillators, 
allowing us a more intuitive interpretation of the differ- 
ent regimes. Thus the changes from recurrent to dormant 
states of tumor cells in the periodical dosage regime can 
be interpreted as a lock of the unstable oscillations of the 
nonlinear oscillator imposed by the linear one. We can 
understand the complex behavior of coupled oscillators 
by representing its dynamics as a function of control pa- 
rameters V and (3 Jl7]Jl8|]. In order to depict it we plot, 
in the parameter space (Vvs/3), those points for which 
tumor growth is uncontrolable. 

In the case - < 1 for a < — , for effective value of doses 
and frecuencies higher than certain threshold, the sys- 
tem can revert from uncontrolable growth to a treatment 
controled population. Namely, for every set of parameter 
values a,k and a, there are threshold values for [3 and V 
which split the parameter space (see Fig 6) into two zones 
corresponding to uncontrolable and controlable growth of 
malignant cells. The Fig 6 was generated for a fixed set of 
initial conditions. Although the behavior of the param- 
eter space for different initial conditions is qualitatively 
the same, the threshold values show strongh dependence 
of the initial conditions. 

From this result and taking into account the meaning of 
the parameters p and V, we can infer that treatment is 
specific for each patient and kind of tumor since thresh- 
old doses values depend on the immunological response of 
each individual, on the malignant cells population at the 
begining of the treatment and also on the rate of prolif- 
eration of the tissue. Besides, the existence of threshold 
values reflects the fact that reaching controlable popula- 
tions of malignant cells is only possible by mantaining a 
minimal dose above certain threshold given by the con- 
tinuous line depicted in Fig 6, which can be well fitted 
by a hiperbolic function. 

For effective doses and frecuencies below these thresh- 
old values, the system behaves qualitatively the same as 
without treatment. However, contrary to this statement 
it is also found that for low frecuencies the growth of ma- 
lignant cells can be controlled in spite of being below the 
threshold values (see, for details, Fig 6). 
There exists a "paradoxical" phenomenon observed in 
experiment and the clinic, consisting in the fact that the 
enhacement of the inmune system with immunotherapy 
stimulates tumor growth [ fl9[ , which could be explained, 
"in some way", by this result, i.e., why, for fixed doses 
burden , the growth becomes uncontrolable at given fre- 
cuencies above those localized in the region of controlable 
growth of malignant cells (Fig. 6). Such values would rep- 
resent an optimal treatment as it reduces doses burden 
and treatment frecuency. 

Also, for higher frecuencies with small V, growth of can- 
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cer cells can be controlled as shown in Fig 6. The obser- 
vance of these optimal values would be important because 
of the negative effects produced when cytokines concen- 
tration reaches above a critical concentration [^0|,^T) . . 
Now, increasing the amplitud of the effective dose value 
for a fixed frecuency above the threshold, a malignant 
cells population reduction is obtained, in spite of an un- 
controlable growth being observed for some higher doses 
burden (see Fig 6). 

On the other hand, setting V to some value and varying 
the frecuency from zero to higher values, different behav- 
iors are reproduced. At zero treatment frecuency, tumor 
cells population is controlable with the lowest values of 
doses burden (Fig. 6). However, for frecuencies different 
from zero, we find zones of recurrent and dormant growth 
of tumor cells. The population of cancer cells controled 
under treatment presents an oscillating behavior (see Fig. 
7) 1- 

For values of the parameters satisfying — < a < 1 (that 
we interpret as a dormant state), malignant cells popu- 
lation can be reduced by increasing effective doses. On 
the other hand, varying the frecuency for fixed effective 
doses values, an oscillating behavior for the population 
of malignant cells is obtained, as in the previous case. 
In all these cases, regrowth of malignant cells takes place 
after treatment interruption ||. This can be easily un- 
derstood if we take into account that the population of 
malignant cells with zero value (x — 0) represents, in the 
mechanical analogue (Eq.[l4]), an unstable point (a po- 
tential maximum, as that shown in Fig. Ib). This means 
that any variation will lead the system towards a minimal 
potential position. Therefore for a residual population 
slightly greater than zero, a regrowth of tumor cells will 
take place after the treatment. 

Analysing the behavior when — > 1 we arrive at the fol- 
lowing results. The range of values with physical sense 
for a, i.e. I < — < a and I < a < — allows only two 

'a a J 

critical points: a stable node and a saddle point. In this 
case dynamics in the phase portrait is the same as that 
without treatment. There are no possible changes in the 
dynamics, so treatment is useless. 

Hence, for — > 1, only initial conditions determine the 
final outcome of tumor evolution, irrespective of the ap- 
plied treatment. 

VI. CONCLUSIONS 

In this work we intend to give a new focus to the 
dynamics of tumor growth in a periodical regime of im- 
munotherapy. We explain such dynamics considering this 
system as two coupled oscillators, namely, the compe- 
tence between the immune system and malignant cells 
analized as a nonlinear oscillator coupled with a linear 
one that simulates the treatment. 

This simple model allowed us to describe all possible 
states in wich a tumour can be found. It also presented 



some of the features found in tumor dynamics, outlined 
by some authors, such as the existence of short term os- 
cillations of tumor size as well as the long term tumor 
relapse. On the other hand, this model gives the depen- 
dence of tumor growth on some parameter values related 
to the treatment: the frecuency and amount of applied 
doses. We conclude from this study that the evolution of 
tumor submitted to immunotherapy has a strong depen- 
dence on these parameter values. 

In some cases, growth of malignant cells can be reverted 
with inmunotherapy treatment. Corresponding thresh- 
old values are obtained for treatment frequency and dose 
above which growth is stopped and malignant cells pop- 
ulation reduced. 

Also, for those inmunological parameters for which a sta- 
ble population of malignant cells exists, the size of the 
dormant tumor can be reduced by increasing dose bur- 
den. 

It was shown, as well, that for certain relation among 
the parameter values, tumor presents a recurrent behav- 
ior with or without treatment. 

In all these cases, when a reduction of the tumor is pos- 
sible, best results are obtained for low constant values of 
the dose. Nevertheless in all cases previously analyzed, 
after the interruption of the treatment, tumor regrowth 
is observed. 

This would confirm the fact that treatment with im- 
munotherapy using cytokines alone is not succesful 
enough in the treatment against cancer Q. Therefore 
another kind of therapy would be required. 
As a way to improve the model, we propose the intro- 
duction of other terms taking into account effects pro- 
duced by stochastic perturbations due to enviromcntal 
conditions ^,^,^1). Some authors atribuit to these per- 
turbations the main cause of possible jumps from stable 
to unstable behavior in tumor growth dynamics. This 
will be considered in a future work. 
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FIG. 1. Potential barriers in which the particle moves, (a) 
a = 1, k = 1.5 and a = 3. (b) a = 2, k = 0.2 and a = 0.25 .. 

FIG. 2. States and its stability for the second fixed point 
as a function of parameters a, k, a. 



FIG. 3. Evolution of malignant cells on time without treat- 
ment, (a) a = 2, k = 0.2 and a = 0.05 with xo = 2.1 and 
yo = 2.7. (b) a = 2, k = 0.2 and a = 0.25 with x = 5.3 and 
yo = 6.7. 



FIG. 4. Phase portraits (lymphocytes population versus 
malignant cells), (a) a = 2, k = 0.2 and a = 0.09. (b) 
a = 2, k = 0.2 and a = 0.25. 

FIG. 5. Saddle point in the phase portrait for values of 
parameters a = 1, k = 1.5, a = 3. 

FIG. 6. Growth behavior of malignant cells with treat- 
ment for a = 2, k — 0.2 and a = 0.05 with 
xo = 5.3 and yo = 6.7 depicted in the parameter space 
(V, 0). Uncontrolable growth (gray points). Controlablc 
growth (white points). Black solid line (hiperbolic function 
V = 0.10478 + 0.00044/(0.05343 + /3) 2 ' 7313 ). 

FIG. 7. Limit cycles for V — 0.25 and different values of 
parameter (3 in the phase portrait. 
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